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We develop the Floquet scattering theory for quantum mechanical pumping in mesoscopic con- 
ductors. The nonequilibrium distribution function, the dc charge and heat currents are investigated 
at arbitrary pumping amplitude and frequency. For mesoscopic samples with discrete spectrum we 
predict a sign reversal of the pumped current when the pump frequency is equal to the level spacing 
in the sample. This effect allows to measure the phase of the transmission coefficient through the 
mesoscopic sample. We discuss the necessary symmetry conditions (both spatial and temporal) for 
pumping. 
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I. INTRODUCTION 



Quantum charge pumpiny3~0 is presently of consider- 
able interest. An experiment by Switkes et al. ^ demon- 
strated that a phase coherent mesoscopic system sub- 
jected to a cyclic two parameter perturbation can pro- 
duce a directed current. 

Coherent quantum pumping is a consequence of the in- 
terference of energetically different traversal paths made 
possible by an oscillating scatterer. The ratio of the oscil- 
lation frequency uj to the inverse time taken for carriers 
to traverserihe-sample r^T^i^efines the operational regime 
of a pumpES^EJ. BrouwerB gave an elegant formulation 
of adiabatic (w ^ ^t^) quantum pumping that is based 
on the scattering matrix approach to low frequency ac 
transport in phase coherent mesoscopic systems devel- 
oped by Biittiker, Thomas, and PretrelS. This approach 
leads naturally to apgeogietrical description of adiabatic 
quantum pumpingpilj 113. The theory predicts that the 
charge pumped during a cycle depends on the area en- 
closed by the path in the scattering matrix parameter 
space. A less formal but more physical picture of an 
adiabatic quantum pump appeals to both quanturn-^Ba- 
chanical interference and photon assisted transport .tSlH. 
The same processes i-ajpe, important for a nonadiabatic 
{lo ^ T^^) pumpEirEalla. These discussions emphasize 
the energetics of the carrier traversal process. 

It is the purpose of this work to develop a theory that 
permits the description of both adiabatic and nonadia- 
batic regimes on the same footing and allows a simple 
physical interpretation. To this end we extend the ap- 
proach of Ref. ^ to the case of large frequencies and 
large pumping, amplitudes. We apply the Floquet scat- 
tering theoryE3~Ej which deals with the scattering matrix 
dependent on two energies (incident and outgoing) . This 
approach leads to expressions for the quantities of in- 



terest in terms of the side bands of particles exiting the 
pump. The side bandsO correspond to particles which 
have gained or lost one or several modulation quanta fiu). 
This approach is complementary to discussion.s based on 
the scattering matrix dependent on two timesOouJ 

The paper is organized as follows. In Sec.|| the gen- 
eral approach to the kinetics of quantum pumps based 
on the Floquet scattering theory is presented. In Sec. Ill 
we apply the general results to the adiabatic case. In 
Sec.^ we calculate the Floquet scattering matrix for a 
particular model - an oscillating double barrier potential 
- and present the results of numerical calculations of the 
pumped charge and the heat currents in both adiabatic 
and nonadiabatic regimes. We conclude in Secj^. 



II. GENERAL APPROACH 

We consider scatterin^lS'l of an incoming flow of 
electrons with energy E dX a. scatterer that oscillates in 
time with frequency uj. |— . 

During the interaction with the oscillating scattererta 
electrons can gain or lose energy quanta fiuj. Hence the 
outgoing state is characterized by the set of energies i?„, 
n = 0,±l,±2,.... 



En — E - 



ihu. 



(1) 



This is a Floquet state. 

According to the Floquet theorem the energy ladder 
Eq. (|l]) gives the fuUpStit of possible energies for outgoing 
particles (see e.g.,EBEj). Thus to describe scattering due 
to an oscillating scatterer we can use the Floquet scat- 
tering matrix Sp- The matrix element SF,af3{Em E) is 
the quantum mechanical amplitude for an electron with 
energy E entering the scatterer through lead p to leave 
the scatterer through lead a having absorbed {n > 0) or 
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emitted (n < 0) energy quanta |n|?iu;. The Greek let- 
ters a, /3, number the leads connecting the sample to Nr 
reservoirs. 

We remark that the negative values En < correspond 
to bound states near the oscillating scatterer. These 
states influence scattering into the propagating {En > 0) 
states but they do not directly contribute to the current. 

Current conservation implies that the submatrix 5^^' 
of the Floquet scattering matrix (corresponding to prop- 
agating modes only) is a unitary matrix 



c{out) 



Mp)f Mp) _ c(p)Q(p)t _ f 



(2) 



In particular, if a current with flux 1 and energy E enters 
the scatterer through lead /3 then current conservation 
implies 



E \SFMEn,E)f = 1. 



(3) 



E„>0 



Another useful condition follows from the fact that if all 
incoming propagating (£'„ > 0) channels are full then 
each outgoing channel has also to be full: 



E E \SFME,En) 
13 E„>0 



1. 



(4) 



Note that usually the Floquet energy E is determined 
within the interval Q < E < fiu). However, for our prob- 
lem, it is convenient not to reduce the discrete set of En 
to this interval and to keep E as the actual energy of 
incident (or outgoing) particles. 

Because of Eq.(|^) we can express the annihilation op- 
erator b for outgoing particles in the lead a in terms of 
annihilation operators a fprJLticoming particles in leads 
/3 = 1, 2, . . . , as followsE3'El 

L(i;) = E E (5) 

13 E„>a 

The operators «„(£') obey the following anticommuta- 
tion relations 

[aiiE),afi{E')] ^ dc.„d{E ~ E'). 

Using Eqs.(|2|) and (||) we see that the operators ba{E) 
obey the same relations. 

Note that above expressions correspond to single 
(transverse) channel leads and spinless electrons. For the 
case of many-channel leads each lead index (a, f3, etc.) 
includes a transverse channel index and any repeating 
lead index implies implicitly a summation over all the 
transverse channels in the lead. Similarly an electron 
spin can be taken into account. 

Now we calculate the distribution function /^°"*^(i?)Q 
= < h]^{E)ha{E) > for electrons leaving the scatterer 
through the lead a. Here < ... > means quantum- 
statistical averaging. Taking into account Eq.(|^) we ob- 
tain 



{E)^Y: E \SF.ME,En)[' 

13 -B„>0 

-^ft\En). 



(6) 



Here /^*"^(i?„) — < a^ij{E)ap{E) > is the distribution 
function for electrons entering the scatterer through lead 



A. Directed charge currents 

Using the distribution function /q""*-* (E) for outgoing 

particles and /q*"'' (E) for incoming ones we can find the. 
directed current /„ in the lead a far from the scatterercS 



Io. = l f dE [fi°-'\E) - fi"^\E)} 



(7) 



The current directed from the scatterer towards the reser- 
voir is positive by definition. Substituting Eq.(^ into 
the above equation, using Eq.(|4|), and making the shift 
E ^ E — nfav we find 



oo 

Tl^dEY: E \SFMEn,E)\'' 
f3 £;„>o 



[ff''\E)-f^r\En)) 



(8) 



Here Eb >o means a sum over those n (positive and 
negative) for which £"„ — E + nfiu > 0. 

Another useful representation for the directed current 
can be obtained if we use Eq.(||) and make the shift 
E ^ E - nfiuj in /i°"*^(£;) in Eq.(0). As a result we 
obtain 



l3^aE„>0 

{\SFMEn,E)\^f;-\E) 



(9) 
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.MEn,E)\^ft\E)} 



IFtoui this expression for the directed current it follows 
that only transmission a ^ (3 (not reflection a = f3) con- 
tributes to the current. In addition Eq.(^) can help us to 
consider the effect of time reversal symmetry (TRS) on 
the pumped current. 

On the one hand, the time reversal t —t (TR) in- 
terchanges incoming and outgoing channels 

[SF,al3{En,Eyf^^^ = S F,f3a{E , En) ■ 

Hence if the TRS is present then Eq.(||) reads 
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(\SF,a0{En, E)\'^ — \SF,ap{E, En) 



(10) 



In the above equation (in accordance with the usual 
pump setup) we suppose that incoming electrons in all 
the channels are described by the same Fermi distribu- 
tion function with temperature T and electrochemical 
potential ^: 

fi"'HE) ^ foiE) ^ 
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1 E-n 



Generally, the Floquet scattering matrix elements for 
transmission with incident energy E to En is not equal 
to the transmission from En to E, 

SF,af3{En, E) ^ S F,al3{E , En) ■ 

^From this we can conclude that even a pump with TRS 
can generate a directed current. If these two scattering 
amplitudes are not equal, there exists the possibility of 
empty states deep below the Fermi-surfaceoO. In this 
case interaction and inelastic effectsM can be expected to 
be especially important. 

However, if the scattering matrix is energy indepen- 
dent on the scale of the order of huj 

SF,al3{En, E) W SF,ap{E), 

then the scattcrcr with TRS can not produce a dc cur- 
rent. The last circumstance is especially important for 
the adiabatic case ^ 0, since the adiabatic scattering 
matrix Eq.(|l^) satisfies the above condition. 

On the other hand, in general, the nonstationary case 
is without TRS and the oscillating scatterer can, in prin- 
ciple, generate a current. In particular, if two parame- 
ters of the scatterer oscillate in time with the same fre- 
quency Lo but with phase lag At/? then the time reversal 
implies the reversal of the sign of A(y9. Thus for a pump 
with A(/5 ^ the scattering problem is without TRS (the 
time reversal symmetry is dynamically broken) and such 
a scatterer can produce a dc current. This effect was 
used far generating a dc current in adiabatic quantum 
pumpsEJ. 



B. Directed heat currents 



By analogy with Eq.(0) we find the directed heatxjp. 
rent lE,a flowing in lead a away from the scattererEjEJ'Ej 
(we suppose that all the reservoirs are at the same elec- 
trochemical potential /i) 



Ie.o. = j^J dE{E - n) 



x{/(™*)(£;)-/^)(£;)} 



(11) 



With the distribution function for outgoing particles 
j{out) Eq.(l) we get 



Ie,c. ^UdEE E (En-fJ.) 



13 B„>0 



:\SFMEn.E)\''[f'^;-\E)~f^^-\En)) 



(12) 



Note that if all the reservoirs are at the same macroscopic 
conditions (electrochemical potential, temperature, etc.) 
then the heat fiow lE,a (at any lead a — IJL . . . , Nr) is 
directed from the scatterer to the reservoircS. That dif- 
fers strongly from the charge current la given by Eq.(|^) 
which, if it exists, can be directed either from the reser- 
voir to the scatterer (at some lead) or vice versa (at an- 
other lead). 



III. ADIABATIC APPROXIMATION 



In this section we use the above formalism to investi- 
gate the limit of adiabatic scattering. 

The general physical notion of adiabaticity applied 
to the scattering problem of interest here is as follows: 
Let us suppose that the time independent problem is 
described by the scattering matrix Si^{E, Xi, . . . ^ Xn) 
which depends on the energy E of incident electrons and 
a set of parameters X^, i = 1,2, . . . , N. Next assume 
that the parameters Xi vary in time: Xi = Xi(t). Then 
we can say that the nonstationary scattering problem is 
adiabatic if scattering of particles incident with energy 
E can be described via the scattering matrix 5*0 with 
time-dependent parameters Xi(t) 



SadiE,t) ^ SoiE,Xiit),. . . ,XNit)). 



(13) 



This approximation is adequate if the scattering matrix 
changes only a little while an electron interacts with the 
scatterer. In other words, the characteristic time scale 
for the change of parameters Xi is much larger than the 
traversal time tt- 

An analogous criterion can be formulated concerning 
the energy dependence of the scattering matrix. To this 
end we consider the adiabatic problem when the param- 
eters change periodically in time: Xi{t) — Xi{t + 2tt /uj). 
In this case we can expand the adiabatic scattering ma- 
trix into the Fourier series 



Sad{E, t) 



E 



Sa^niE)e' 



(14) 



where 
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xSo{E,X,{t),...,XN{t)). 



(15) 



The Fourier harmonics 5o,n define the amphtudes of side 
bands for particles traversing the adiabatically oscillat- 
ing scatterer with initial incident energy E. Thus we 
can construct the adiabatic Floquet scattering matrix as 
follows 



SF,ad{En, E) — SF,ad{E,E^n) — Sa^n{E). 



(16) 



The adiabatic Floquet scattering matrix consists thus 
for each incident energy of a block of dimension of 
{2NrnmaxY' where Nj. is the number of leads and rimax 
is the maximum number of side bands needed for an ac- 
curate description of the quantities of interest. 

The above equation allows us to formulate the follow- 
ing adiabaticity criterion:. If the Floquet scattering ma- 
trix Sp changes only a little when the energy E changes 
by Umax^Lj then the adiabatic approximation can be ap- 
plied and Sf ~ SF,ad- Note that because of Eq. (p^) the 
same criterion can be applied to the scattering matrix 
So{E). 

Let us now calculate the directed charge current lad 
and a heat flow lE,ad in the adiabatic limit. Substitut- 
ing Eq.(|l6|) into Eq.(^ we flnd the current flowing in 
lead a under an adiabatic change of parameters (we put 
/("^ = /o): 



rJdEi- 



dE 



(17) 



xEE"I^O,o/3,n(i?)|^ 

13 n 



We see that the current lad is due to photon-assisted pro- 
cesses (the current depends on the intensity of side bands 
which is proportional to |S'o,Q/3,n(£^)P)- Therefore even 
in the "adiabatic" limit pumping is, strictly speaking, a 
nonadiabatic phenomenon. 

Note that formally Eq. ( p^ ) is obtained for huo <^ ksT. 
However it gives the correct answer (in the adiabatic 
limit) for an arbitrary ratio of the frequency lu to the 
temperature T. 

Using Eq.(p^ we can rewrite Eq.(p^) as follows 



27r/(.i 



Iad,a ^ij^ J dtJdE 



dfojE) 
dE 



(18) 



(^SliE,t)) 



This is Brouwer's formulao. We emphasize that Eq.(|l7|) 
and Eq.(|l|) express a (nonadiabatic in nature) pumped 
current in terms of the scattering matrix So for a " frozen" 
scatterer with time-dependent parameters (see Eq.(|3|)). 
This is correct if the scattering matrix 5*0 can be taken 



to be energy independent on the scale of the order of 
Timax^- The last conclusion is correct irrespective of 
the amplitude of the oscillating parameters. But Umax 
might very well depend on whether one is in the weak 
amplitude limit or in the strong amplitude limit. 

Next consider the heat flow lE.ad,a produced by the 
adiabatically oscillating scatterer in the lead a. Sub- 
stituting the adiabatic scattering matrix Eq.(|l6|). into 
Eq.(|l2|) we find 



i / rfeEE|5o,o/3,„(M + e)P 

— oo /3 n 

(,+ li|-)si„h(5^) 



(19) 



cosh(.^) 



+cosh ( 



At zero temperature (or if the temperature T is less than 
the relevant energy scale for the scattering matrix Sq) the 
above equation can be greatly simplifled. In this case the 
heat flow can be expressed in terms of an on-shell, time- 
dependent adiabatic scattering matrix 5*0 (/i,t): 



2tt/u 



lE,ad,a — ^ J dt 



x(^^^^(A*,^)) 



(20) 



This expression coincides with that given by Avron et al. 

m 

In the adiabatic limit it is possible to express the quan- 
tities of interest in terms of the stationary scattering ma- 
trix. This permits the very general and useful expres- 
sions discussed in this paragraph. In contrast, in the 
non-adiabatic case, the information contained in the sta- 
tionary scattering matrix is by definition not sufficient. 
As a consequence, in order to address the non-adiabatic 
case, we have now to consider a very specific example. 



IV. CHARGE AND HEAT FLOWS PRODUCED 
BY AN OSCILLATING DOUBLE BARRIER 



To apply the general formalism of Sec.|| beyond the 
adiabatic approximation we investigate a simple model |^. 
It is a one-dimensional scatterer consisting of two delta 
function barriers oscillating with frequency ld and located 
at X = —L/2 and x = L/2. 

To calculate the Floquet matrix we consider scatter- 
ing of electrons with energy E coming from x = —oo. 
The time-dependent Schrodinger equation for an electron 
wave function \I'(x,i) reads 
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(21) 

V{x, t) = Vi{t)S{x + f ) + V2{mx - I), 
V^{t) = Va, + 2Vu cos{ujt + ip^), i ^ 1,2. 

Appendix A gives the exact solution of this modeL We 
now use this solution to calculate the charge current flow- 
ing through (as well as the heat current flo wing from) the 
oscillating double barrier potential (see Eq.(|2^)) connect- 
ing two reservoirs with the same temperature T = and 
electrochemical potential /i. We use the units 2m — h — 
e = 1. To be definite we consider the current Ii {Ie,i) 
flowing to the left [x — > — oo). The currents are calcu- 
lated, using Eq.(||) and Eq.(|l2|) with the Floquet scat- 
tering matrix determined as shown in Appendix A. We 
compare these currents with the adiabatic currents Iad,i 
and lE,ad.i using the Brouwer formula Eq. (|l8|) and its 
analog Eq.(|2^). 

For reference we write down the scattering matrix 
So{E) which we need to calculate the adiabatic currents: 



Sq — 



Here k — 

(1 - A)e-^'=^; 




2^ sm{kL) 



1 



(22) 



C + 2f sin(A:L) 



V,m/h' (j 



1) 



1,2); C = 

P1+P2 QjT 



^(e V . ^ — 

interest is a comparison of the exact Floquet scattering 
matrix and adiabatic theory in the limit of small and 
large frequencies. Furthermore, it is interesting to com- 
pare the symmetry conditions for pumping for the two 
theories. 



A. Adiabatic limit: 



^From our calculations it follows that in the limit of 
w — > the adiabatic current lad gives a good approxima- 
tion for the pumped current irrespective of the amplitude 
of the oscillating potentials Vu . However the criterion of 
adiabaticity depends strongly on the ratio Vu/fuv. 

As we have seen, the adiabatic approximation is valid 
if the scattering matrix Sq is energy independent on the 
scale of the order of nmaa^i'-^- Heren,„oa; (specified below) 
is the number of side bands with noticeable amplitude ex- 
cited by the oscillating scatterer. Let us denote by SE the 
energy scale over which the scattering matrix 5*0 changes 
significantly. For example, close to a resonance SE is of 
the order of the width 5 of a resonance. Then the adia- 
batic approximation is valid if rimaxfi'jJ ^ SE. Hence the 
larger the number Umax of excited side bands the smaller 
the frequencies for which the adiabatic approximation is 
valid. 
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FIG. 1. Relative change of the pumped current 
[lad — lii-^)]/ lad as a function of the pump frequency u) for 
three values of an oscillating potential Vu = V12 ~ 0.02 (solid 
line); 1 (dotted line); and 2 (dash-dotted line) close to the 
transmission resonance. The frequency is measured in units 
of the width of the resonance 5 = 0.339. The parameters are: 
L = 2-k; fi = 24.167; Voi = V02 = 20; (pi - f2 = tv/2. We use 
the units: 2m — h — e — 1. 

As it is well knownil the number rimax of excited side 
bands depends on the ratio of the amplitude of the os- 
cillating potential Vu to the frequency lu. In the small 
amplitude limit Vu — > only the first side bands are ex- 
cited: Umax = 1- But for stroug pumping the number 
iT-rnax is large and the adiabatic approximation is valid 
only at smaller frequencies. As a consequence for a given 
finite frequency to the deviation of the actual pumped 
current /(w) from the adiabatic one lad increases with 
increasing pumping amplitude. This fact is illustrated in 
Fig.0. 

Before concluding this subsection we would like to em- 
phasize the following. For strong pumping the adiabatic 
approximation is still valid at sufficiently small frequen- 
cies despite the excitation of a large number of side bands. 
The Fourier harmonics of the scattering matrix Sq define 
the amplitudes of side bands for strong (Eq.(|l^)) as well 
as for weakE3 adiabatic pumping. Thus the scattering 
matrix Sq completely defines the kinetics and, in par- 
ticular, the quantum statistical correlation properties of 
an adiabatic pump. The current noise of a weak (small 
amplitude) pump was considered in Ref. |2^. The noise 
in a strong amplitude adiabatic pump will be presented 
elsewhere. 



B. Large frequency limit 



At large frequencies u! > 5 (where 5 is the width of 
a resonance) the pumped current / differs considerably 
from the adiabatic one lad- 
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FIG. 2. The adiabatically pumped charge 5Qad ~ 
(in units of an elementary charge e) during a period 
(left panel) and the adiabatic heat current lE,ad (in units 
of huj'^/{ATv)) (right panel) as a function of the Fermi en- 
ergy /i. The parameters are: L = 2007r; Vbi = V02 = 20; 
Vii = V12 = 10; <pi - <fi2 = 7r/2. 



It is known that for a double-barrier structure at w — s- 
the pumned dc current /aduEj and the generated heat 
flow Is.ac^ show a resonance-like behavior as a function 
of the Fermi energy n (see Fig.^) . The pumped current 
and heat peak when the Fermi energy is close to the 
transmission resonance. This is because in the adiabatic 
limit only the particles close to the Fermi level contribute 
to the charge (Eq.dl^)) and energy (Eq.(po|)) transfer. 



With increasing pumped frequency 3> 5 the particles 
within a wider energy interval come into play. The depen- 
dence on the Fermi energy is smoothed away. However a 
resonance-like dependence on the frequency uj arises. 



In Fig.|| we depict the dependence of the charge 6Q 
= {2tt)I/lu pumped during a cycle on the frequency co. 
The pumped charge peaks when the energy quantum huj 
equals one (or several) level spacings A of a double bar- 
rier structure. In the example used for numerical cal- 
culations the level spacing A near the Fermi energy is 
constant with good accuracy. The smaller peaks corre- 
spond to many [n — 2,3, . . .) photon processes for which 
nhoj = mA (m = 1,2,. . . ). We see that the single photon 
processes dominate over the many photon processes. 



FIG. 3. The charge 5Q — 2ttI/uj pumped during a cycle 
as a function of the frequency lu (in units of the separation 
A = 0.0417 between the transmission resonances). The Fermi 
energy is ^ = 17.423. The parameters are the same as in 
Fig.2. 

At arbitrary {nhui ^ mA) but large {huj ^ 6) fre- 
quencies, due to interference (inside the double barrier 
structure) only the main component {n ~ 0) of the Flo- 
quet state Eq.(23) has a significant amplitude between 
the barriers. This component corresponds to the eigen- 
function of the time independent problem (with barriers 
Vol and Vb2)- The side bands (n ^ 0) do not participate 
in the transmission through the system (more precisely, 
their contribution is small) . As a result the pumped cur- 
rent is greatly reduced compared to the adiabatic case. 
It should be noted that in the adiabatic case the side 
bands do contribute to the transmission since for w ^ 5 
they lie at the same transmission resonance as the main 
component. 

On the other hand at some particular values of a fre- 
quency fiio = toA the substates '0±fc {k = 1,2,3,.. .) of 
the Floquet state Eq.(23) become large and additional 
channels for the transmission through the system open 
up. This leads to the increase of the pumped charge (see 
main peaks in Figj^). 

Interestingly, as shown in Fig.^), the pumped current 
reverses sign as a function of frequency. These sign rever- 
sals can be understood in the following way: If two po- 
tential barriers oscillate then the pump effect arises as-a 
consequence of an interference between two amplitudest^. 
The first amplitude Ai corresponds to particles which 
propagate through the double barrier structure and ab- 
sorb (emit) the energy quantum fiLD (or several quanta 
nhoj) at the first barrier. The second amplitude A2 cor- 
responds to the same propagation with absomtion (emis- 
sion) at the second barrier. As it is knownldfl the adia- 
batically pumped current is an odd function of the phase 
difference of the oscillating potentials A(/7 — ^pi - f2 (see 
Eq. (pTl) ) . In fact this is the phase difference of the two 
amplitudes Ai and A2- In the nonadiabatic case there 
exists an additional contribution coming from the spatial 
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phase difference 



A(^a; ~ {ki - k2)L. 



Tlius at large frequencies the pumped current depends 
on the sum Aip + Aip^- 

To clarify the appearance of A(px let us consider par- 
ticles with energy E going from the left to the right and 
absorbing the quantum Hlo. The amplitude Ai corre- 
sponds to particles which absorb an energy tko close to 
the left barrier (Vi) and traverse the system at energy 
E + hu>. Hence fci ~ \/E + Twj. On the other hand the 
amplitude A2 corresponds to particles which absorb an 
energy Uu) close to the right barrier (V2) and thus tra- 
verse the system at energy E. They thus propagate with 
wave vector k2 ^ \fE. Because k\ 7^ ki there is a phase 
difference A^p^ between the amplitudes A\ and Ai ■ Note 
that in the adiabatic case (i.e., at w — > 0) Atp^ vanishes. 

Now we show that in the case under consideration the 
phase difference Akp^ is at the origin of the sign reversal 
of the pumped current at consecutive peaks (see Fig.||). 
Close to the m*'* main peak of /(w) the frequency is deter- 
mined by fiLO ~ mA. Taking into account that only the 
particles with energy close to the resonance contribute to 
the transmission and mainly single photon processes are 
important we can estimate A<^x as follows. If the am- 
plitude A\ corresponds to the propagation through some 
resonance (say Ei+m) then A2 corresponds to the propa- 
gation through the resonance Ei. At the transmission 
resonance it is fc^'^L = In the additional phase difference 
is Aifx = rmr. Thus the sum Aip + Aipx changes by tt as 
we pass from one peak ta another (Am = 1) thus giving 
rise to the sign reversalc3. 
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FIG. 4. The heat current Ie (in units of huj'^ / (An)) as a 
function of the frequency uj. The parameters are the same as 
in Fig.3. 

In addition interference inside the double barrier struc- 
ture manifests itself as a resonance-like dependence of the 
heat production (see Fig.g). We can see that Ie{'~ij) con- 
tains additional peaks at hto = (m + i)A. This fact 



emphasizes the striking difference between the processes 
which are pesponsible for charge pumping and for heat 
productionEa. 



C. Quantum pumping and broken symmetry 



Our results show that the symmetry of the scattering 
problem is important for the pump effect. 

The quantum pump effect arises due to the different, 
interfetiag excitation histories of carriers traversing the 
sampleE3. To "extract" directed currents from the uni- 
form environment the left-right symmetry (LRS) for car- 
riers traversing the sample needs to be broken. There 
are, at least, two ways to break the LRS. The first one 
is breaking the spatial symmetry (SS). The second one is 
breaking the time reversal symmetry (TRS). 

We restrict our considerations to systems which in the 
absence of the time-dependent perturbations needed for 
pumping are in an equilibrium state. In particular this 
means that without the presence of magnetic fields (etc.) 
the system is time reversal invariant. Note that the time 
independent scattering problem is insensitive to the pres- 
ence (absence) of the spatial symmetry of the scatterer 
(in the sense that the transmission probability is invari- 
ant under the spatial inversion x — > ~x irrespectively of 
the SS of the scattering potential). 

Interestingly, the symmetry conditions for pumping de- 
pend on the frequency. In the adiabatic limit w — > 
(see Sec. [II) the scattering problem is fully characterized 
by the scattering matrix 5*0. In contrast to the conduc- 
tance, adiabatic pumping is sensitive to the symmetry 
of the scattering matrix. To have an adiabatic pumping 
effect we need to break both the spatial symmetry and 
the time-reversal symmetry. 

It is important to distinguish the symmetry of the 
equilibrium problem and the symmetry of the full prob- 
lem. In the adiabatic case, the symmetry of the time- 
independent problem is irrelevant: What matters, is that 
the system in presence of perturbations breaks both the 
spatial symmetry and the time reversal symmetry. In 
many examples (for instance the two barrier problem con- 
sidered above) it is not possible to break the time-reversal 
symmetry without at the same time breaking the spatial 
symmetry. However, examples which are spatially sym- 
metric and have broken time-reversal invariance can be 
constructed. For instance, we could consider the problem 
of two barriers oscillating in synchronism and as a sec- 
ond perturbation which oscillates with a phase lag choose 
the potential between the two barriers. This would be an 
example of a two parameter problem which does not gen- 
erate a pumped current. The important point we would 
like to emphasize is that the symmetry conditions for the 
non-adiabatic pump are different. 

In the strongly nonadiabatic limit the scattering 
problem is described by the Floquet scattering matrix 
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5'^(^;(o«*)^i;(i»)) dependent on both the incident E^^"^ 
and outgoing energies. In general, the Floquet 

scattering matrix is sensitive to the spatial symmetry of 
the scatterer. If the scattering problem is without SS 
(i.e., when the scattering potential is not invariant under 
the change x —x) then 

and the oscillating scatterer can pump a dc current (see 
Eq.(^). Thus to obtain a pump effect in the strongly 
nonadiabatic limit there are, at least, two possibilities. 
Either the spatial symmetry or the time reversal symme- 
try (or both SS and TRS together) have to be broken. 
The role of spatial symmetry breaking for the nonadia- 
batic pump effect was emphasized in Ref. |lq. 

The above reasoning is illustrated in Fig.H where the 
dependence of a pumped current on the frequency is 
shown for three generic situations. First, two potentials 
oscillate with the same amplitude but out of phase (solid 
line). In this case the perturbation breaks both the spa- 
tial symmetry and the time reversal symmetry. Second, 
only one potential oscillates (dotted line). The TRS is 
present but the SS is dynamically broken. The system 
is able to pump a dc current like in the previous case. 
Third, two potentials oscillate in phase and with the same 
amplitudes (dash-dotted line). Now both the TRS and 
SS are preserved and the system under consideration does 
not exhibit a directed current. 




-0.1 

2.5 2.7 2.9 3.1 3.3 3.5 

CO/A 

FIG. 5. The charge 5Q pumped during a cycle as a func- 
tion of the frequency uj. Three cases are presented: (i) two 
potentials Vii — V12 = 10 oscillate out of phase ifi - ip2 ~ 
7r/2 (solid line); (ii) only single potential oscillates: Vu — 10, 
V12 = (dotted line); (iii) two potentials Vu = V12 ~ 10 os- 
cillate in phase (pi — ip2 (dash-dotted line). The parameters 
are the same as in Fig. 3. 



V. DISCUSSION AND CONCLUSION 



In this work we have developed the Floquet scattering 
matrix approach to parametric pumping in phase coher- 
ent mesoscopic systems of noninteracting electrons. We 
have calculated the distribution function, the dc current 
and the heat flow at arbitrary pumping amplitudes and 
frequency. 

The Floquet scattering matrix describes naturally the 
existence of the side bands En = E + nhu of parti- 
cles leaving the pump. These side bands correspond to 
nonequilibrium particles generated by the pump which 
carry the heat from the oscillating scatterer to the reser- 
voirs and (under spme conditions) transfer charge be- 
tween the reservoirg23. If the pumping amplitude is small 
only the first (n = ±1) side bands are excited. But for 
strong pumping the number of excited side bands is large: 
nmax S> 1. This number together with the frequency de- 
fines the energy scale UmaxJi'-^ characteristic of the pump 
problem. In particular this is important for the analysis 
of the conditions under which the adiabatic approach to 
pumping is valid. If the scattering matrix is energy in- 
dependent on the scale of the order of rimax^i^ then the 
adiabatic approximation can be applied. In this case the 
elements of the Floquet scattering matrix are given by 
the corresponding Fourier components of the stationary 
scattering matrix Sq (see Eq.(|lq)) with parameters taken 
to be dependent on time. 

The existence of the pump effect is directly related to 
the symmetry of the scattering problem. In the adia- 
batic case (uj ~* 0) only a scatterer without spatial and 
time reversal symmetry can produce a directed current. 
This conclusion applies irrespectively of the (stationary) 
spatial symmetry of the scatterer. On the other hand 
in the nonadiabatic case (at large pumping frequency) 
to achieve pumping we need a scatterer with either bro- 
ken spatial or time reversal symmetry. Hence at large 
operating frequencies even the scatterer with a single os- 
cillating parameter can show a pump effect if only the 
spatial symmetry is brokerJla. 

To emphasize the main physics underlying pumping 
we have considered an exactly solvable model. Namely 
a one-dimensional scatterer consisting of two oscillating 
delta-function barriers Eq.([2l|) separated by a distance 
L. 

If two barriers oscillate out of phase (broken TRS) the 
basic process leading to the pump effect is an interference 
of two quantum mechanical amplitudes. These ampli- 
tudes Ai and A2 corresponds to particles which traverse 
the scatterer in the same direction but gain (or lose) a 
modulation energy in the vicinitjz. either of the first or of 
the second barrier, respectivelytj. This becomes clearer 
at large frequencies w because of the following. When the 
energy quantum huj equals one (or several, say m) level 
spacing A of a double barrier structure the pumped cur- 
rent peaks (see Fig.^. The particles with energy close to 
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the resonance of a double barrier structure give the main 
contribution to the pump current. Thus the amphtudes 
Ai and A2 correspond to traversal of the system through 
different (say, {I + mf^ and l*^) resonances. This gives 
rise to an additional phase difference l^ipx = (fci — k2)L 
between the amplitudes Ai and A2- For the resonant 
double barrier considered here the phase difference be- 
tween consecutive resonances is tt and thus lS.ifx = mn. 
As a result the consecutive peaks of the pumped current 
have opposite sign (see Fig.g). 

We conclude that if the time reversal symmetry is bro- 
ken (i.e., Aip 0), the nonadiabatic pumped current 
contains direct information on the phase of the trans- 
mission coefficient through the mesoscopic sample. We 
emphasize that here we have a single-connected geometry 
in contrast with the case where the mesoscopic sarogle (a 
quantum dot) is embedded in one arm of the ringcZI. 

We have also considered the effect of a spatial asym- 
metry. If only one potential (say, Vi) oscillates then the 
spatial symmetry is (dynamically) broken. In this case 
the (nonadiabatic) pump effect arises because the proba- 
bihty for exciting side bands |S'f,q/3(-E„, E)]"^ depends on 
the direction (i.e., on the order of indexes a and /3) (see 
Eq.(|)). The particles going from the left to the right 
first absorb (emit) the energy nhui — En — E and then 
pass through the double barrier structure (at energy En). 
Thus the corresponding probability is 

\SFMEn,E)\^ ^T{En). 

Here T{E) is a transmission probability. On the other 
hand the particles going from the right to the left pass 
the system before they absorb (or emit) the energy in the 
vicinity of the oscillating barrier I. Hence we have 

\SFS2{En.E)\^ ^T{E). 

Thus the pumped current Eq.(^ depends on the differ- 
ences T{E) - T{E ± nhuj). 

Note that for a system with equidistant spectrum (that 
is the case under consideration) the transmission prob- 
ability T{E) is periodic in energy E with the period of 
A. Hence if Hui = mA then it is T{E) « T{E ± fiuj) and 
the current is close to zero in accordance with Fig.|^ (a 
dotted curve; m = 3). 

We have presented a scattering theory of quantum 
pumping based on the Floquet theorem. This approach 
allows a description of the kinetics of both adiabatic and 
nonadiabatic quantum pumps. In particular the quan- 
tum statistical correlation properties (noise) can be con- 
sidered in analogy with the weak amplitude adiabatic 
quantum pumpcJ. The investigation of noise is especially 
important in tba., context of quantized (i.e., noiseless) 
charge pumpinglj-llfrLJ and will be presented elsewhere. 
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VI. APPENDIX A: OSCILLATING DOUBLE 
BARRIER. AN EXACT SOLUTION 



In this Appendix we determine the Floquet scattering 
matrix of the Schrodinger equation (|2^) of two oscillat- 
ing delta function potentials located at a; = — i/2 and 
X — L/2. To find the Floquet matrix we consider scat- 
tering of electrons with energy E coming from x = —fo.. 
To solve ( ^T|) we use the Floquet functions methoda. 
Since the Hamiltonian H depends on time, the system 
has no stationary eigenstates. However, the Floquet the- 
orem tells us that because the Hamiltonian is periodic in 
time the eigenstates of Eq.(^T|) can be represented as a 
superposition of wave functions with energies shifted by 
nhuj: 



*is(a;,i) 



-iEt/h 



(23) 



n— — oc 



Away from the points x = —L/2 and x — L/2 the func- 
tions 'ipnix) are a superposition of plane waves 



(24) 



where 



2m 

kn = d (E + nhcu) 

with Re[kn] > and /m[fc„] > 0. At a: 
X — L/2 we have the boundary conditions 



-L/2 and 



0,i) = *B(-|+0,t), 



*i5(f -0,i) = *B(f +0,t). 



d-fE{x,t) 
dx 



dx 



(25) 



a'S'E{x,t) 

dx 



dx 



x = 4f-0 



x^'^+O 

In the case under consideration (i.e., when an electron 
with energy E comes from x = —00) the functions ipnix) 
are 



tpnix < 0) = SnflS'^'^'' + Rn 



-ikjiX 



4>niO <x<L)^ A„e*'="- + B„e-*'="^ (26) 
V'„(x > L) =r„e^'="^ 
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Here the coefBcients Rn and T„ for propagating modes 
(for which £"„ > 0) are the ampUtudes of reflection from 
or transmission through the double barrier system ab- 
sorbing [n > 0) or emitting [n < 0) an energy |n|?iijj, 
respectively. 

Substituting Eq.(p^) and Eq.(p6|) into the boundary 
conditions Eq.d23) we obtain 



(27) 



Rn — (A-n — Snfi) 6 ' " + i3„ 



T — A -I- R p-ik^L 



where the coefficients and i3„ are subject to recursive 
equations. It is convenient to represent these equations 
in a matrix form 



An 
Bn 



(28) 



An+l 



An-1 



Here we have introduced the matrices 



{ikn -poi)e' 



Un 



-P02e 



{ikn-po2)e 2 



(29) 



Dn+l = I 1 , (30) 



Piie 




where the parameters are: pij = Vijm/h , i = 
J = 1,2. 

To solve the equation (^^ we have generalized the 
method used in Ref. ^ for a single oscillating delta- 
function potential. 

First of all we consider Eq.(p8|) at n > 0. We suppose 
that there exist matrices X„ such that 



An 
Bn 



Xr, 



An- 
Bn- 



(32) 



Then substituting the above constraint into Eq.(p8|) we 
obtain a simple recursive equation for the matrices X„ 
(n > 0) 



Xn = {U„ 



In- 



(33) 



Using these matrices we can express all the coefficients 
An, Bn {n > 0) in terms of Aq and Bq only 



An 
Bn 



Ao 
Bo 



^X„...Xi[ ] , n>0. (34) 
Further we consider n < and introduce the matrices 



An 
Br, 



An+l 
Bn+l 



(35) 



The corresponding recursive equation for the matrices Yn 
is 



Y, 



D 



n+1- 



(36) 



Using these matrices we express the coefficients An, B„ 
(n < 0) in terms of and Bq: 



An 
Bn 



Yn 



^0 

Bo 



n<0. 



(37) 



As a final step we consider Eq.([2^) at n = 0. After simple 
manipulations we find the coefficients Aq and Bq: 



Ao 
Bo 



Uo - DiXi - I^iY^i 
ikoe~^^°' 



(38) 







Thus, using the solutions of the recurrent equations 
Eq.(p3[) and Eq.(|3^) we can calculate all the coefficients 
An, Bn (see Eq.(^), Eq.(|^, and Eq.(||)). Note that 
in each particular case we need to take into account only 
the limited number |n| < Umax of side bands and thus 
we can put Ar„^^^+i w and Y^nmai-i ~ 0- 

The coefficients i?„ and r„ can be calculated with 
the help of Eq. (^ . For the propagating modes {En = 
E+nTiui > 0) of interest here, these coefficients determine 
the elements of the Floquet scattering matrix 



\SF,2liEn,E)\'^ — J^\Tn\'^- 



(39) 



Here the indexes 1 and 2 correspond to the left and right 
reservoirs, respectively. To obtain the matrix elements 
Sf,22 and S'_f,i2 we need to solve the same problem with 
plane waves coming from the right. 
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